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I. Introduction 

Owing to non-linearity of the Schwinger-Dyson (S-D) equation^"'' of the fermion self- 
energy, together with running of the gauge coupling constant being included in, it is 
scarcely possible to obtain an analytic solution of the equation. As a result, one usually has 
to solve it by numerical method^~^°. However, one could still get some analytic solutions 
of the S-D equation if some linearization approximation to the equation is made. The 
resulting analytic solutions will be quite useful for discussions of chiral symmetry breaking. 
In this paper, we will derive and expound the general exact analytic solutions of the S-D 
equation of fermion self-energy in the linearization approximation. 

Consider a theory with vectorial gauge and chirally-invariant effective four-fermion 
interactions. When running of the gauge coupling constant is taken into account, the S-D 
equation of the fermion self-energy T,{jP') in the Landau gauge, after Wick rotation and 

angular integration, will have the following form^^: 
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(1.1) 

mo is the bare fermion mass, A is the ultraviolet (UV) momentum cut-off, C2{R) is 
the eigenvalue of the squared Casimir operator of the "color" gauge group in the rep- 
resentation R of the fermion field ip with the dimension d{R), h is the strength of the 

chirally-invariant four-fermion interaction — ^'^^ ^^(max(p2, A;^)) is the 

conventional approximation of the running gauge coupling constant g'^{(p — kY) defined 
by 

5'^(p^) if > k"^ 
f{k^) if A;2>p2 
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Eq.(l.l) will be reduced to that 
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In Eq.(1.4) we have used a continuous Ansatz of the running gauge coupling constant 



where 



and 



-g\q') = l//3or(g^) 



r(g^)=ln(^ + 



Po 



(1.5) 
(1.6) 

(1.7) 



llC2(G)-^4T(%)iV; 
/ 

with the standard denotations in gauge theory. In the flavor sum Z!/ in Eq.(1.7), besides 
the fermions corresponding to ip, all the lighter colored fermions will also be included in. 
The scale parameter /i is optional and the parameter ^ is required to be greater than 1 so 
as to avoid the infrared (IR) singularity of g^{q^)- 

The integral equation (1.4) is equivalent to the following differential equation 



u;{x)J:"{x) + [u;'{x) + l]E'(x) 



E(x) 



t{x) [x + E2(a;)] 



(1.8) 
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together with the IR boundary condition 



E'(0) = 



2(lnOS(0) 



;i.9) 



and the UV boundary condition 



1 + —xt{x] 



uj{x)E'{x) + E{x) \ = mo(A) 

:A2 



where 
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(1.10) 
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We emphasize the following points: 1) The non-linearity of Eq.(1.8) is embodied in the 
term in the right-handed side of the equation. 2) The IR boundary condition (1.9) comes 
from the x — > limit of the equation 



E'(x) 




y + E^{y) 



(1.13) 



which is a result of Eq.(1.4). Since Eq.(1.13) is valid for any x, we can define the IR 
boundary condition by Eq.(1.13) at some non-zero value of x, instead of at x = 0, if 
E'(a;) is calculable at that value of x. This fact will have a close bearing on the actual 
solution of the following linearized S-D equation. 3) The term in Eq.(1.4) relevant to the 
four-fermion interactions contains no x hence appears only in the UV boundary condition 
(1.10) rather than in Eq.(1.8) itself. 

In Sect. II we will state the linearization approximation for Eq.(1.8) and a necessary 
change of the IR boundary condition of the solution. In Sect. Ill the derivation of the 
general analytic solutions of the linearized S-D equation will be given in detail and in 
Sect. IV we will conclude with some discussions on forms and features of the independent 

ll. I'lie linearization approximation 

The non-linear S-D equation (1.4) or (1.8) has no analytic solution. Some analytic 
solutions could be obtained merely in the linearization approximation of the equation. To 
see how to make this approximation, we first consider the IR and UV asymptotic solutions 
of Eq.(1.8). 

In the region where x is small, Eq.(1.8) becomes 



Suppose E(a;) has the form 

E(x) = x^{ao + Qix + a2X^ + ■■■), when x is small 
and substituting Eq.(2.2) into the IR boundary condition (1.9) i.e. 

limE'(x)E(x) = --^ 
x^o ^ ^ ^ ^ 21n^ 



(2.1) 



(2.2) 



(2.3) 
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we obtain that 

alsx"^"'^ + (2s + l)aoaia;^*' + \{s + l)al + 2{s + l)aoa2 



lim I c 



X 



2s+l 



21n{ 
(2.4) 

The only possibihty to satisfy Eq.(2.4) with a real ao is to set s = 0. With this result, 
substituting Eq.(2.2) into Eq.(2.1) we will have 



2ai+6a2X+12a3X +• • •+ 



ao 



1 

ao 



+ ai] X + 



„ai a, \ n 
+ 3^ + ^ - aa + 
at, ao I 







(2.5) 

Thus we may express ai, a2, as • • • by means of ao = E(0) and write down the solution 
of Eq.(2.1) in small x region which is consistent with the IR boundary condition (1.9) as 
follows: 

b h f . h 



E(x) = E(0) <^ 1 



-X 



2(lnOS2(0) 6(lnOS4(0) 
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(2.6) 



12(lnOS6(0) 

where E(0) is a finite constant. 

In the region where x oo, we may have two possible assumptions: 1) S^(a;) > x 
and 2) S^(x) < x. If E^(x) > x is assumed, the asymptotic form of Eq.(1.8) will become 



that 
X ( 1 



1 



Inx 



E"(x) 



Inx 



E'{x) 



b' 



(Inx)E(x) 



, b' 



1 + A 



with A < 1 (2.7) 



where the ratio a;/E^(x) at x — * oo in the denominator of the right-handed side of Eq.(2.7) 
has been replaced approximately by the constant A. This replacement will not change 
the essential behavior of the asympototic solution. Eq.(2.7) remains to be a non-linear 
equation. The exact form of its asymptotic solution is in general unknown. However, we 
can always suppose a typical form of the solution, such as the form of conventional power 
function, so as to examine whether the assumption that E^(a;) > x when a; — > oo could 
be consistent with Eq.(2.7) or not. Thus we set that 

J:{x) ^C{lnxYx' (2.8) 

where C is a constant, then Eq.(2.7) may be reduced to 

C {-r(r - l)(\nxY-^ + r(r - 2s - l)(ln x)""' + [r(2s + 1) - s^](\nxy-^ + s(s + l)(lnxY} x''^ 



b' n \-r-l -s 
= - — [Inx) X 



(2.9) 



When s 7^ 0, —1, the leading term in the left-handed side is Cs{s + l){\nxYx^ ^ whose 
equality to the right-handed side requires that 



-r - 1, s - 1 = -s, Cs{s + 1) = -b'/C 



(2.10) 



and they give that r = —1/2, s = 1/2 and = — 46'/3. The results demand E^(a;) ~ 
x/ \nx, contradictory to the assumption E^(a;) > x. When s = or —1, the leading term 

in the left-handed side is C[r{2s + 1) — s^]{\n xY"^ x^"^ whose equality to the right-handed 
side requires that 



r - 1 = -r - 1, s - 1 = -s, C[r(2s + 1) - s^] = -b'/C 



(2.11) 
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and they give that r = 0, s = 1/2 and = 4b'. The result s = 1/2 is obviously opposite 
to the presupposition s = or — 1. Therefore, the assumption E^(a;) > x when x —>■ oo 
can not be consistent with Eq.(2.7), at least this is true for the supposed asymptotic form 
(2.8) of the solution. Alternatively, we may assume that S^(x) < x when x — > oo. In this 
case, Eq.(1.8) approximately becomes that 

^ ^ x^"ix) + (2--^) J:'ix) = -TT^^ (2.12) 



lax J V hixj (Inx) x 

Noting that Eq.(2.12) has now been linearized. We may still substitute the trial solution 
(2.8) into Eq.(2.12) and obtain the algebraic equation 

[r(2s + 1) - s^]{\nxY-^x'-^ + s{s + l)(lna;)V-^ = -b{\nxY-^x'-^ (2.13) 

where the non-leading terms such as (In x)^~^, (In x)'^"^ etc. have been neglected. Eq.(2.13) 
could be satisfied if the term with (lna;)'^a;*~^ is removed by setting 

s(s + l) = (2.14) 

i.e. s = Oors = — 1. Thus we will obtain the UV asymptotic solution of S(a;) with the 
form 

EW^^(.nij)"V^(,ni,)" (2.15) 

The fact that the UV asymptotic form (2.12) of Eq.(1.8) has been linearized indicates that 
the non-linearity of Eq.(1.8) is important only in the IR region. Actually, the linearization 
approximation of Eq.(1.8) is a good one not only in the UV region but also in the mediate 
momentum region where x is not very small. The approximation can be made by replacing 
E^(a;) in the denorminator of the non-linear term in Eq.(1.8) by E^(0) and this will result 
in the following linearized version of Eq.(1.8) 

u;{x)^"{x) + Wix) + mx) = -^J^ll^ (2.16) 



Generally, the E^(0) in the right-handed side of Eq.(2.16) may also be replaced by some 
constant which will be viewed as the average value of S(a;) over small x region. So 
the linearized Eq.(2.16) will be qualitatively correct for small x where non-linearity is 
important. 

To make it be easy to solve Eq.(2.16) analytically we take two further assumptions. 
One of them is to suppose that 

S^(0)=e/^^ (2.17) 

This is permissive since both the IR parameter ^ and the scale parameter fi are undeter- 
minated theoretically. Of couse, the constraint that ^ > 1 will demand that S^(0) > fi^. 
Another assumption is that the coefficient uj{x) [Eq.(1.12)] in Eq.(2.16) is replaced ap- 
proximately by 

u;{x)c^^^^{x + ^t,') (2.18) 



This approximation is valid if 



X 

0^ 



1 + 



> 1 (2.19) 
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The condition (2.19) could impose some constraints on ^. Let x/^ji"^ > r then Eq.(2.19) 
will become that 

The inequahty (2.20) is obviously valid for r > 1 and it may be changed into that 

1 

e < e— for < r < 1 (2.21) 

r + 1 ^ ^ 

Hence, in the case with r < 1, ^ will have an upper bound. The upper bound will decrease 
as r goes down. Although for some values of r we could obtain the constraint on ^ which 
is allowed by phenomenenlonogy (for example, if r = 3/4 then the constraint will be 
1 < ^ <^ 11. 477), in the following we prefer r = 1 to r < 1, i.e. wc will use the assumption 
(2.18) continuously down to the scale x — ^^'^ without the need to consider any upper 
bound of ^. 

Under the assumptions (2.17) and (2.18), the linearized equation (2.16) may be reduced 
to the following form that 

^^^^ {X + i^')T.'\x) + [2- -^^-,] T.\x) = (2.22) 



1 + t(x)^ ^'^ ' \ [1 + t{x)Y\ T{x){x + iiJ?) 

The solution of Eq.(2.22) must submit to the UV boundary condition (1.10) and some IR 
boundary condition as well. Since Eq.(2.22) is inapplicable in the region where x < ^/i^, 
we will not be able directly to use the IR boundary condition (1.9) at x — > 0. Instead, as 
was indicated in Sect. I, we may now use Eq.(1.13) and the asymptotic solution (2.6) of 
E(a;) in small x region to give a new IR boundary condition at x — ^//^, i.e. 



E'(x)L_...,^bi- 



^"^''^ dx \xt{x) 



where we have also made the assumption that the solution (2.6) of S(a;) in small x region 
will be approximately extended to x — ^jji^. The use of the IR boundary condition (2.23) 
instead of Eq.(1.9) will allow us to be able both to solve Eq.(2.22) exactly and analytically 
and to include partially the IR non-linearity of Eq.(2.10) in. 

It is noted that the scales of physical chiral symmetry breaking will generally be in the 
region where x > so it is reasonable to use Eq.(2.22) to discuss such kind of problem. 
The assumptions that the approximation (2.18) and the solution (2.6) arc extended to 
X = ^fi^ respectively from above and from below will at most numerically affect the IR 
boundary condition (2.23) of Eq.(2.22) and not produce an essential impact on physical 

iTf.^^So^utions of the linearized equation 

We will find out the exact indepedent solutions of the linearized equation (2.22). 
Instead of x, let r to be the new variable, then Eq.(2.22) may be changed into that 



:E"(r) + 



1+ ' 



(1 + r)'. 



S'(r) + 6^=0 (3.1) 



1+T 

Assume the exact solution of Eq.(3.1) have the form 

E(t) = e-T7(r) (3.2) 

where s and r arc constants to be determined. Substituting Eq.(3.2) into Eq.(3.1) and 
dividing each term of the equation by e*'^T''+^, we will obtain the differential equation 
satisfied by /(r) as follows: 

(r + 1)/"(t) + [(2s + 1)t + 2(r + s + 1) + 2(r + 1)/t]/'(t) 
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+{s{s + l)r+[2s{r + l) + s^+r + b] + [r{r + l) + 2s{r + l) + 2b]/r+[r{r+l) + b]/r''}f{T) = 

(3.3) 

The UV (r oo) asymptotic solutions of Eq.(3.1) may be obtained from Eq.(3.3) by 

taking /(r) = constant and setting the first and the second term in the coefficient of /(r) 
to be equal to zeroes as well as neglecting all the terms with the orders 1/r and above. 
The results are denoted by 

S|^7eg(T) = and = e-W'-' (3.4) 

Eq.(3.4) is obviously consistent with the UV asymptotic form (2.15) of S(a;). In order to 
find out the exact solutions of Eq.(3.1) with the UV asymptotic forms (3.4), we must set 
the first and the fourth term in the coefficient of /(r) in Eq.(3.3) to be equal to zeroes, 
i.e. 

s(s + l) = (3.5) 

and 

r(r + l) + 6 = (3.6) 

Eq.(3.5) will make S(a;) have the correct UV asymptotic form (3.4) and Eq.(3.6), whose 
left-handed side is just the coefficient of can also be obtained from the small r 

limit of Eq.(3.1) by setting E(r) ~ r*", hence it certainly contains the small r behavier of 
Eq.(3.1). Since the r given by Eq.(3.6) takes complex values, the corresponding solutions 
of S(t) will include complex constants. However, a physical solution may be a real linear 
combination of the two real components of a complex S(r). 

Eq.(3.5) gives that s = and s = —1. We will discuss the two cases respectively. 

1) s=0. In this case, the solution (3.2) becomes E(°)(t) = r7(°)(r). Then from 
Eq.(3.3), the equation of /*^°^(t) can be changed into that 

{z' - z)f'y\z) - [z' - j{z - l)]/(°)'(^) - {az - b)f'\z) = (3.7) 

where 

z = -T (3.8) 

7 = 2(r + l) (3.9) 

a^b + r (3.10) 

Eq.(3.7) could be solved by means of confluent hypergeometric (Kummer) function^^. 
Noting that if \z\ — >• oo, then Eq.(3.7) becomes 

zf'^"{z) + (7 - z)f''y{z) - af'\z) = (3.11) 

which is just the confluent hypergeometric equation with the solution f^^^ = iFi{a; 7; z). 
Hence, the solution of Eq.(3.7) may be written to be that 

f'\z)^aiF,{a;r,z) + 9{z) (3.12) 

where a is a constant. Substituting Eq.(3.12) into Eq.(3.7) and considering that 1^1(0;; 7; z) 
obeys Eq.(3.11) we obtain that 

{z'^ — z)g" — \^ + 7(1 — z^g' + (6 — OLZ^g — az \F\{a.\ 7; z) — a{b — a) iFi{a; 7; z) (3.13) 

It is seen that the left-handed sides of both Eq.(3.13) and Eq.(3.7) have the same form. 
In view of iFi{a; 7; z) has been used in f^^\z) we may assume that 

g{z) = cz iFi'{a;r,z) (3.14) 
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where c is a constant. Substituting Eq.(3.14) into Eq.(3.13), we will meet the terms with 
iFi"{a;j; z) and iFi"{a;'j;z). However, they can be changed into some combinations 
of the terms with only iFi{a;j;z) and iFi{a;j;z) by means of Eq.(3.11). In this way, 
Eq.(3.13) will be reduced to that 

c Q - 1^ ^ iFi{a; 7; z) - ca iFi(q;; 7; z) = az iFi'(q;; 7; z) - a{b - a) iFi{a; 7; z) 



which demands that 



a — 



7 



1 c = 0. 



— {b — a)a + ac = 
The a and c have non-zero solution if and only if 



a 



1 (6 - a) = 



(3.15) 



(3.16) 



which is obviously satisfied based on Eqs.(3.9), (3.10) and (3.6). As a result, we have 
a = | — l = rifc=lis taken. From this we will have the solution of Eq.(3.7) with the 
form 

f^'\z) = (^--l) ,F^{a;r,z) + z^^,F^{a;r,z) (3.17) 
and the corresponding 



Q - 1) iFi(q;; 7; -r) - iFi (« + 1; 7 + 1; -r 



(3.18) 



Depending on the UV asymptotic form, the physically-relevant real solutions will have 
the following two kinds of forms: 



yiO) 
irreg 

regV 



r)^}= 4°)}^~'"^'[(i"0 iFi(a;7;-r)-^riFi(a + l;7 + l;-r) 
where one of two roots of the algebraic equation (3.6) of r 



r = — - — irj, V — ~ 7 



+ c.c. 

(3.19) 

(3.20) 



has been definitely chosen (chosing the other root r = + if) has no difference because 
the physical solutions depend on only the two real components of Ti^^\t)). The complex 
constants Af^ and A^") will be determined respectively by the UV asymptotic conditions 



and 



reg 



reg ' 



By using the asymptotic formula of confluent hypergeometric function 
iFi(a;7;-r) ^ 



12 



r(7)^-r^_^.|a-7^ r(7) 



r a] 



r(7 - a) 



(3.21) 
(3.22) 

(3.23) 
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we can write 



(t) 

^irrcg \ ' / 



) 1 r(7) ^ ■ 
/ r(7 - a) + 



The conditions (3.21) and (3.22) will lead to that 

,(0) • ^ 



Via) 



and 



teinh(27rr/)|r(a)|^r(7) 



^(0) ^ ,r(a) sm(7ra) 



(3.24) 
(3.25) 

(3.26) 



r(7) sinh(27r?7) 

2) s=-l. In this case the solution (3.2) becomes E(~^)(r) = e~^r*'/*^~^)(r) and from 
Eq.(3.3), /^~^-*(t) submits to the equation 

t(t + l)f-^y\r) - [r^ - (7 - 2)t - 7] /("'^'(r) - (ar + 7 - b)f^-^^ (r) = (3.27) 

where 



a — 1 + r — b 

The solution /(^^H'^) 

can be expressed by iFi(q;;7;t). Assume that 
/(-i)(r) = a iFi(a;7;r) + cr iFi'(a;7;r) 



(3.28) 



(3.29) 



Substituting Eq.(3.29) into Eq.(3.27) and considering that iFi(q;;7;t) satisfies the equa- 
tion 

T iFi"(a; 7; r) + (7 - r) iFi'(a; 7; r) - a ^F,{a- 7; r) = (3.30) 

we obtain that 

[-a + {b + a)c]T iFi{a; 7; r) + [{a + b- 7)0 + ac] iFi{a; 7; r) = 
which demands that 



-a + {b + a)c =0 

(a + & — 7)0 + ac — 

The condition on which a and c have non-zero solution 

a-F (6-Fa)(6-Fa -7) = 



(3.31) 



(3.32) 



is certainly satisfied in view of Eqs.(3.28),(3.6) and (3.9). Hence if taking c—1 then we 
will have a — b + a — r+l = 7/2 and 



/(-i)(r) = I iFi(a;7;r) iFi(«;7;r) 



and correspondinly, 



^ iFi(a; 7; r) + -r iFi{a + 1; 7 + 1; r) 
2 7 



(3.33) 



(3.34) 



The real solutions nl^^lUr) and -'^^(t) can be expressed by 



irreg V / 



r_-2-"7 



e T 2 



^ iFi(a; 7; r) + -r + 1; 7 + 1; r) 

2 7 



+ C.C. (3.35) 



where the complex constants and will be determined respectively by the UV 

asymptotic conditions 



and 

By using the asymptotic formula of confluent hypergeometric function 

r(7) 



iFi(q;;7;t) 



r(7) 

Fa 



r(7 - a) 



(3.36) 
(3.37) 

(3.38) 



we get from Eqs.(3.36) and (3.37) that 



^(-1) ^^^_par(Q;) sin(7rQ;) 



and 



r(7) sinh(27r77) 

TT r(Q:) 



6sinh(27rry)|r(a)|2r(7) 

IV. Some discussions of solutions 



(3.39) 



(3.40) 



We will show that among the four real solutions S-^^g and T^'f^^ (s = 0,-1) derived 
above, only two ones are linearly independent, consistent with the general conclusion that 
a second order differential equation has two independent solutions. 
In fact, E(°)(t) and E(~^)(t) differ by only a constant. Denote 



EW(r) =r7^°^(r) 



then Z'^^) will obey Eq.(3.7) or 



T(r + l)/(°)"(r)+ r^ + 7(r + l) /^"^'(^) + + = 



Now write 



/(o)(r) = e--/i(r) 



then it is easy to verify that hij) obey the equation 

t(t + l)h"{T) - fr^ - (7 - 2)t - 7I h'{T) - (ar + 7 - 6)/i(t) = 



(4.1) 

(4.2) 
(4.3) 

(4.4) 



which has the same form as Eq. (3.27) of /(-I) (r). It follows from this that /i(r) oc /(-^)(r). 
Hence 

E(°) (t) = e-W'h{T) oc e- VV^-^) (r) = E^"^) (r) (4.5) 

i.e. Tj^^\t) and E'^~^)(t) differ at most by a constant and both are linearly dependent. 

This conclusion may also be reached directly from the explicit expressions of E('^)(t) 
and E(-i)(t). By means of the property of confluent hypergeometric function 



iFi (a; 7; -r) = e iFi (7 - a; 7; r) 



(4.6) 
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we can express E(°^(t) in Eq.(3.18) by 



r) = e V 



(1 \ ^ 
- - Ij iFi(7 - a; 7; r) - -r 1^1(7 - a; 7 + 1; r) 

Prom Eqs.(3.9), (3.10) and (3.28) we see that 

7 — a = q; + 1 
By Eq.(4.8) and the recurrence relation^^ 

7 iFi{a + 1; 7; r) = T iFi{a + 1; 7 + 1; r) + 7 iFi{a; 7; r) 
we may rewrite Eq.(4.7) in the form that 



(4.7) 

(4.8) 
(4.9) 



E(o)(r) 



e-r^r (^^ _ 1^ ^. ^) + i (1 _ _a_^ ^ ^^^(^ + 1; 7 + 1; t) 

>i(a; 7; r) + f f T iFi(a + 1; 7 + 1; r) 



7 



(4.10) 



iFi(a;7;T) + fr iFi(a + 1; 7 + 1; r) 



'E(-i)(t) 



where Eqs.(3.6), (3. 9), (3. 10), (3.28) and the expression (3.34) for ^^(r) have been used. 
Eq.(4.10) shows that E*^°)(r) and E'^~^)(r) differs by only the constant (7 — 2)/7 indeed. 
By comparing Eq.(4.7) with Eq.(4.10) and using the relation (4.6) we see that the equality 



^ \ a 'y — 2 

^ - 1 iFi(a; 7; -r) r ^F^{a+l■, 7+I; -r) = ^ e"^ 

z / 7 7 



^ iFi(a; 7; r) + -r iFi(a + 1; 7 + 1; r) 
2 7 

(4.11) 



is valid. From this equality it may be verified that the real solutions in Eqs.(3.19) and 
(3.25) satisfy the relations 



sfciiWand E(°)(r) = EL-^)(r) 



(4.12) 



It is emphasized that in the verification of Eq.(4.12), the equahty (4.11) hence the rela- 
tion (4.6) plays a key role. Since the constants Af'\ Al,^\ A\^^^ and A^"^-* are determined 
through the UV asymptotic form of TjI^^^^{t) and S*!^^(r) (s = 0,-1) this means that 

the relation (4.6) must keep to be valid also in the UV asymptotic region. It is just this 
requirement that determines the selections of the phase factors in the asymptotic expan- 
sions (3.23) and (3.38). 

In summary, the linearized Schwinger- Dyson equation (3.1) of the fermion self-energy 
in the theory with running gauge coupling constant has two linearly independent real 
solutions which can be considered as the the complex solution T,^^\t) or 'E^~^\t) or the 
pairs of real solutions (Ei„eg(r) = e£^(t) = eH^(t), Ercg(r) = ^^%{t) = E[-^\t)). 
The physical solutions should be some linear combination of Eirreg(T) and Ereg(T) with 
the combination coefficients to being determined by the IR and UV boundary conditions 
(2.23) and (1.10). The application of these physical solutions to problem of chiral sym- 
metry breaking will be discussed elsewhere. 
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